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ABSTRACT 


This  report  describes  the  development  of  a method  for  numerically  computing  the 
inviscid  transonic  flow  field  over  an  airfoil  with  a leading-edge  slat  or  a trailing-edge  flap, 
the  approach  is  to  solve  the  full  inviscid  irrotational  flow  equations  about  two-element  air- 
loil  systems.  The  methodology  consists  of  the  development  of  a suitable  computational 
i)lane  and  grid  system  through  the  application  of  a series  of  conformal  mappings.  The 
appropriate  set  of  equations  and  boundary  conditions  are  derived  in  terms  of  a smoothly 
vai'ying,  single-valued  reduced  potential  function  through  analytic  removal  of  all  singularities 
in  the  computational  domain.  A stable  and  accurate  relaxation  procedure  is  established  for 
the  numerical  solution  of  the  governing  equations.  The  method  is  applied  to  a variety  of 
transonic  supercritical  two-element  airfoil  configurations.  Results  are  presented  depicting 
•he  surface  pressure  distribution,  streamline  patterns  in  the  physical  and  computational  do- 
mains, and  Mach  number  contours. 
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1.  I\TK(Jl)l  C 'l  lo:- 


111  ol’dci  for  an  aircraft  to  maneuver  effeetiwlv  in  liu  ri  an.  oriie  .-|iee(l  raiim  . it.'  iwn;; 
must  generate  high  lift  coefficients  without  incurring  exci'.'.'i  = i-  drag  or  Imlleting.  1 he  re- 
cent development  of  supercritical  wings  can  cnalile  ;i  de-Mgiu-i  to  meet  this  speeiiicaiioii. 
However,  this  requirement  will  often  degrade  the  airerali’.-  iierlormance  at  cruising  siieeds 
with  larger  than  optimal  drag  coefficients.  The  implementation  of  liigh  lilt  devices  at  tran- 
sonic speeds  offers  the  possibility  of  greatly  enhancing  tlic  maneuvering  capabilities  of  mod- 
ern aircraft  without  compromising  their  cruising  efficiency.  This  possibility  has  been 
proven  in  the  last  lew  years  by  the  installation  of  slats  on  the  F-4  and  the  jiositive  test  oi  a 
slatted  wing  on  the  F-14  aircraft.  The  performance  of  thee  aircraft  in  managing  climbs 
and  turns  at  transonic  speeds  was  remarkalily  improved  bv  the  presence  of  the  slats,  even 
tliough  these  configurations  have  not  been  shown  to  be  optimal  ly  any  means. 

The  aerodynamic  designer  currently  lacks  an  analytical  tool  to  design,  or  even  analyze, 
transonic  airfoils  with  high  lift  devices.  Furthermore,  the  paucity  of  experimental  data  cur- 
rently maizes  it  difficult  to  determine  what  can  be  achieved  with  these  maneuvering  devices. 
Also,  the  e.xijeri mental  gathering  of  data  on  such  coidigu rations  would  lie  tremendously  ex- 
pensive in  light  of  the  number  of  configtirations  that  would  need  to  be  tested  and  the  high 
speeds  and  Reynolds  numliers  required  in  a wind  tunnei  lest.  A theoretical  tool  for  the  anal- 
ysis of  the  transonic  flow  over  a two-element  airfoil  system  would  lie  a first  valuable  step  in 
aiding  the  designer  in  his  task  liy  cutting  down  on  the  number  of  configurations  to  be  tested 
and  providing  insight  into  the  flow  phenomena  that  are  present  at  high  speeds. 

4'his  report  describes  the  development  of  a method  lor  numerically  computing  the  in\  is- 
cid  transonic  flow  over  an  airfoil  witli  a leading-edge  slat  or  a trailing -edge  fla]).  In  general, 
these  flow  flchls  are  difficult  to  obtain  analytically  because  of  the  eomjilicated  geometry  of 
the  multiply  connected  domain,  bmall  disturbance  approximations  (such  as  tliose  used  for 
this  problem  in  Ref.  1)  do  not  appear  to  lie  adequate  since  the  interaction  of  the  flow  between 
the  airfoils  will  provide  large  perturbations  to  the  flow  field. 

In  recent  years  the  a])plication  of  mLxed-flow  relaxation  techniques,  introduced  l>y  Mur- 
man  and  Cole  (Rel.  2),  has  made  possilile  the  numerical  computation  of  inviscid  transonic 
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flows  over  a variety  of  geometrical  shapes  in  lx)th  two  and  three  dimensions.  These  methods 
are  generally  based  on  the  assumption  of  irrotational  flow  and  solve  either  the  full  potential 
equation  or  an  appropriate  form  of  the  small  disturbance  equation.  For  two-dimensional 
flows  in  particular,  accurate  and  efficient  solutions  to  the  full  potential  equation  have  been 
ot)tained  for  transonic  flows  over  airfoil  sections  (Refs.  3,  4),  over  axisymmetric  bodies  (Refs. 
5,  (J),  and  over  nacelles  (Refs.  7,  8). 

In  this  report  these  relaxation  techniques  are  applied  to  determine  the  flow  about  an  air- 
foil with  a slat  or  a flap  at  transonic  speeds.  The  approach  is  to  solve  the  full  inviscid,  ir- 
rotational flow  equations  about  two-element  airfoil  systems.  The  methodologj'  consists  of 
the  1)  development  of  a suitable  computational  plane  and  grid  system,  2)  evaluation  of  an  ap- 
propriate set  of  governing  inviscid  eciuations  and  boundary  conditions  in  terms  of  smoothly 
varying,  single-valued  functions  in  the  computational  domain,  and  3)  establishment  of  a stalRe 
and  accurate  numerical  ijrocedure  for  the  solution  of  the  governing  equations. 

An  abbreviated  version  of  this  analysis  was  presented  at  an  International  Symposium 
(Ref.  9)  and  this  report  is  designed  to  provide  the  details  of  the  method.  Arlinger  (Ref.  10) 
has  recently  developed  independently  a similar  inviscid  analysis  of  this  problem. 

'I'he  method  will  provide  solid  groundwoi-k  for  the  development  of  a computational  tool 
that  includes  v'iscous  corrections,  and  one  that  actually  designs  two-element  airfoil  systems: 
a numerical  computation  that  generates  airfoil  ordinates  when  a desired  pressure  disti'ibu- 
tion  is  specified  over  paiT  or  the  whole  of  the  configuration. 
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2.  MAPPINGS 

A crucial  step  in  the  clevelupmenl  oi  a I'inile  diiterciu  c ir.cih'xl  i.i  .impute  flows  ovei' 
complicated  geometries  is  to  develoii  a suital)le  gi’id  system.  It  is  lu^;lily  dcsiralile  to  b.avc 
the  geometric  contours  aligned  with  a coordinate  lino  in  oriier  to  a\oid  iinerjiolations  and  ex- 
trapolations in  appljing  the  surface  boundary  conditions.  It  is  aiso  convenient  foi'  e.xternal 
flow  problems  to  map  the  infinite  physical  domain  to  a finite  computational  space  in  onier  to 
apply  accurately  the  far-field  boundary  conditions.  I'urthermori',  the  maiipinn^  shouki  con- 
centrate grid  lines  in  regions  of  steej)  flow  gradients  such  as  airfoil  Icaiiina  aiui  trailing 
edges  and  in  the  slot  formed  IxUween  the  main  airfoil  *'.nd  the  slat  or  flap. 


In  our  approach,  we  use  analytic  and  numerical  conformal  maiiping  methods  to  trans- 
form the  infinite  domain  external  to  the  two-element  aiifoil  system  to  the  annular  region  bc- 
bveen  two  concentric  circular  rings.  The  outer  ring  corresponds  to  tlic  main  airloil  surtaie 
and  the  inner  ring  to  the  secondary  airfoil  surface  tflap  or  slat).  Infinity  in  'he  physical 
plane  is  mapped  to  a single  point  with  the  circular  annulus  in  the  compuiatiom.  1 domain.  The 
mapping  method  follows  from  the  work  of  Ives  Utef.  11)  and  utilizes  a seciucnce  of  fixe  con- 
formal transformations,  three  analytic  and  two  numerical.  We  briellx’  summarize  IjcTow  the 
details  of  the  mapping  as  they  hax  e been  used  in  our  a]jproach. 


The  first  step  consists  of  mai)ping  the  main  airfoil  to  a near  ciiclo  using  a von  Karman- 
Trefftz  transformation,  as  illustrated  in  I'igs.  la  and  lb.  Denoting  the  physical  plane  as 
Zj  - xi  +iyi,  the  first  mapped  plane  Z2  becomes 

D/kj 


Z2  - s _ / zi  - ziT  y 
Z2  r ^ ' Zl  — ZlN  ' 


(1) 


where  Zif  is  the  locatioi  of  the  trailing  etlge,  Zm  is  located  inside  the  airfoil  at  a point  mid- 
way lictween  the  nose  and  its  center  ol  curxature,  (,'i  = 2 — Tj/t  with  t,  being  the  inchuk'd 
tr.'xiling  edge  angle  of  the  main  ;urfoiI  and  (Zir  — ^ :n  )/-'■!  • 


'I'hc  second  step  is  a Thcoilorsmi-Garrick  tixmsfonnaiion  utilizing  fast  Touricr  trans- 
forms to  map  the  main  aiifoil  nctir  circle  to  an  exai't  circh  . 1 he  lU'Xt  maiiiied  plane.  Z . 

illustrated  in  Tig.  Jc  is 


( A,  iDJZ^ 
f \ 


/ - /. 


A)  Z,  PLANE 


MAIN  AIRFOIL 


B)  Zp  PLANE 


Fig.  1 Sequence  of  Conformal  Mappings 


C)  Z3  PLANE 


MAIN 

AIRFOIL 
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cificle 


D)  Z4  PLANE 


MAIN 

AIRFOIL 


.AT 


Fig  1 Sequence  of  Conformal  Mappings  (continued) 


where  Z2c  is  the  location  of  the  centroid  of  the  near  circular  ima^e  of  the  main  airfoil  in  the 
Z2  plane  and  Aj  and  Bj  are  constants  wliich  are  determined  iterativel}-  as  discussed  in  Ref.  11. 


The  thii'd  step  of  the  mappiii”'  is  an  analytic  transformation  of  the  secondary  airfoil  to  a 
near  circle,  which  keeps  the  image  of  the  main  airfoil  a circle  (but  of  different  radius).  'I'hc 
transformation  to  the  Z.i  plane  is  written  as 


(Z4  - ZjT)(Zj  - Z4t)  ^ 
(Z4  — Z4N)  (Zl  - Z4N) 


tZo  - Z-^t)  (Z3  - ZiiY 
(Za  - ZoN)  {Zj,  — Z3n)„ 


1 /«2 


Id) 


where 

ZjT  = I/Z2T 
= I/Z3N 
Z4T  = 

24N  = l's/^4N 

K'2  = 2 - T2/7T 

and  i’3  is  the  radius  of  the  circle  corresponding  to  the  main  airfoil  in  the  Z.;  plane  aiul  T2  is 
the  included  trailing  edge  angle  of  the  secondary  airfoil.  I14  a manner  sin^ilar  to  the  Zi  — Z2 
mapping,  the  singular  points  Zyj  and  Z^n  are  placed  at  the  trailing  edge  and  miilwa\-  Ijctween 
the  nose  and  its  center  of  curvature  as  illustrated  in  Tig.  Id.  I'he  unJinown  constants  in  Tcj. 
(3)  are  Z4X,  Z.^  and  r,.  Our  procedure  for  the  evaluation  of  these  constants  ditfers  from 
Ives  (Ref.  11)  and  is  detailed  in  the  appendix. 


The  fourth  step  of  the  mapping  consists  of  using  a ijilincar  transformation  to  |)lace  the 
new  circular  image  of  the  sccondaiw  airfoil  at  the  center  of  main  airfoil  circle  a>  illuslixitei! 
in  Tig.  le.  The  mapping  takes  the  form 


where  the  procedure  for  cominiting  the  constants  A,  B,  and  C is  found  in  Ret.  11. 


The  fifth  transformation  of  the  two  concentric  shapes  to  two  cii’cular  rings  is  pi  idormed 
through  a second  application  of  a Thcodorsen-Garrick  mapihng  of  the  form 


\ ^ 

Z5  = Ze  cxp-iDo  + S 

/ i=0 

= ..oie 


(-  Cj  • il)j)(r3  Z,;)'  • (Cj  )(>V  ZJ'  1 1 . 


(3) 


As  shown  in  Tig.  If,  taking  Z^  = re*  , the  main  airfoil  lies  exterior  to  the  eirele  r = 1,  the 


secondary  airfoil  lies  interior  to  r - r,  and  the  point  of  inlinitv  is  at  r i-.,  ('  0. 


i 
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A further  set  of  transformations  X = X(6)  and  Y = Y (r)  is  used  to  obtain  the  final  com- 
putational domain.  It  is  convenient  to  have  the  point  of  infinity  r = r«  located  at  a fixed  grid 
point  in  the  computational  domain,  since  boundary  values  must  be  applied  there.  Hence,  we 
take  Y (r«)  =1/2  and  define  a parabolic  transformation  by 

r - rg  = (Ai  Y + A2>Y  (G) 

where 

Ai  = 2(1  - 2r..  + r,) 

A2  = 4r«>  — 3rg  — 1 

The  main  airfoil  is  located  at  Y = 1,  the  secondary  airfoil  at  Y = 0,  and  the  point  at  infinity 
is  at  Y = 1/2. 

In  the  azimuthal  direction  we  introduce  a mapping  X = X(6)  in  order  to  give  some  con- 
trol on  the  concentration  of  grid  lines  and  to  locate  the  trailing  edges  at  grid  points  in  the 
computational  space.  This  latter  requirement  will  allow  an  accurate  imposition  of  the  Kutta 
condition  at  the  trailing  edges.  The  mapping  consists  of  two  steps;  first  to  produce  desired 
mesh  spacings 

6 = E sin  0 + F sin  20  + G (7) 

where  the  constants  E,  F,  G are  arbitrax’y,  subjected  to  the  constraints 

E < 4/3  F < 1/6 

The  final  step  of  locating  each  trailing  edge  and  the  point  of  infinity  at  grid  lines  is  obtained 
through 

27t|^0  (0)  - 00  J = X + Bj  sin  27tX  + Bo  (cos  27tX  — 1)  (8) 

where 

(01  — Xi)  sin  7rX2  (^2  - X2)  sin  irXi 

sill  ttXi  sin7r(Xi  — X2)  sin  nXs  sin  tt(X2  — Xi) 

(01  — Xl)  cos  7tX2  ^ (O2  — X2)  cos  7tX2 

sin  ttXi  sin  7r(Xi  — X2)  sin  11X2  sin  ti-(X2  — Xi) 

and  01  = 27t  1^0 (^Ti)  - 0qJ,  02  = (6x2)  “ d’o  = 0(6)  and  Cji  and  612  are  the  trailing 

edge  locations  of  the  main  and  secondary  airfoils,  respectively.  The  values  Xi  and  X2  cor- 
respond to  fixed  mesh  points  for  the  trailing  edges  and  X = 0 corresponds  to  the  location  in 
the  0 plane  (0  = 0)  ol  the  point  of  infinity. 


h 


The  final  computationai  domain  is  sketched  in  Fig.  2.  In  this  plane  a unlturm  gi  id  pro- 
duces the  mesh  distribution  in  the  annular  donuiin  shown  in  Tig,  2 and  in  llie  pliysical  domain 
shown  in  Figs,  4 and  u tor  typical  flap  and  slat  configurations.  The  mapping  produces  a 
grid  where  each  airfoil  surface  is  a coordinate  line,  the  trtiiling  edges  oecmr  on  mesh  points 
and  with  a lugh  density  of  grid  lines  in  the  slot  region  and  near  all  stagnation  points.  Al- 
though the  mapping  procedure  is  quite  complicated,  our  computer  program  to  calculate  the 
coefficients  of  all  the  terms  generally  requires  less  thuij  10  seconds  on  an  IHAl  dTO/lGb. 


SECONDARY  AIRFOIL 


Fig,  2 Computational  Domain 
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Fig.  3 Annular  Domain 


3.  MATHK  MATIC^\L  FORMULA  TION 


The  governing  equations  tor  the  mviscid,  irrotational,  compressible  flow  at)out  a two- 
element  airfoil  s3’stem  are  written  in  the  annular  domain  7.^  = re^*  using  the  metric  11  ol  the 
conformal  mappings,  Lqs.  (1)  — (5),  in  terms  ol  a potential  function  4?,  as 


(a-  - v=)4>,r  - ~ -tTe  " ~ 

■ (a“  • u^)^4>r  • (v^  • u^)^vllr  = 


1 (V  - 1) 


(1  - u"  - v") 


where 


V - — <i’ 

' 1!  '■ 


u - 


and  M«  is  the  free  stream  Mach  number  and  ■>  the  ratio  of  speciiic  heats.  I he  boundart'  con- 
dition on  the  surface  of  the  airtoils  is 

<l>r  = ti  on  r = I'j  and  r - 1 HI* 

At  the  airfoil  trailing  etiges,  since  the  metric  !)  goes  to  zero,  tlie  Kutta  condition  that  re- 
quires there  be  no  flow  around  the  trailing  edges,  may  be  expressed  as 

= 0 at  both  trailing  edges.  (it!) 

Singularities  in  the  above  formulation  are  seen  to  arise  for  several  reasons.  First,  the 
metric  of  the  mapping  l)ccomes  unl)ounded  at  inlinity  (r  = r„..  6 - u) . It  can  be  shown  tliat 


near  infinity 


d/g  (/.^,  1 ») 


^ as  — r. 


and  K is  a known  complex  constant  which  can  l)i'  taken  to  be  K kjC  The  metric  mat’  then 


be  regularized  by 


H = H 


(14) 


1 

where 

r s I Zg  - r„  I ^ = r"  - 2r^r  cos  6 ri  cos- 6 (15) 

and  II  is  a smooth  bounded  function  which  goes  to  unity  at  Z = r^. 

Next,  the  potential  function  itself  becomes  unbounded  and  multivalued  near  the  [joint  ol 
infinity  in  the  computational  domain.  One  contriljution  to  the  singular  potential  comes  Irom 
the  behavior  of  the  uniform  flow  near  infinit}',  which  can  Ijc  shown  to  be  ol  the  form,  for 

<l>j(r,  9)  = Real  \ [ = - !^  [r  eos(6  a — ^2)  ~ <''os(a  - ko)]  (lb) 

f Zg  — r„  ) 1 

where  a is  the  angle  of  attack  of  the  free  stream  velocity  vector. 


A second  contribution  to  the  singrdar  potential  comes  from  the  multivalued  nature  ol  tlie 
circulatory  flow  near  infinity  in  the  annular  domain.  In  taking  a closed  circuit  aljout  - I'x, 
the  potential  must  jump  by  2tt  times  the  circulation  about  each  airfoil.  The  solution  lor  the 
circulatory  flow  potential  valid  near  infinity  is  found  as  a solution  to  the  Prandtl-Glauert 
equations  (see  Ref.  12)  and  is  transformed  to  the  computation  domain  as 

4>.(r,  ^)  = - (Tj  Tj)  tan‘»  [>' 1 - tan  p]  (17) 

where 


= ko  — O'  ■ 7: 


tan 


/ r sin  9 \ 

Vr  cos  9 — 1’  / 


(18) 


and  r,  and  T2  are  equal  to  2n  times  the  circulation  constants  about  the  main  and  secondary 
airfoils,  res[)ectively. 


To  obtain  a single-valued  reduced  potential,  another  term  T3  must  Ijc  added  to  4>2  so  that 
any  closed  contours  alxjut  individual  airfoils  will  produce  the  required  circulation  jump. 

This  is  obtained  through  a lei  11 

-1>3  = - F20  (19) 

A reduced  potential  function  G(r.  9)  may  now  be  defined  which  remains  bounded  and 
single-valued  throughout  the  entire  annular  domain  as 

(i(r,  6)  = 4>  — <l>j  — 'ho  — 4>3 
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(20) 


Tlie  governing  equation  lor  the  reduced  potential  is  then  obtained  by  introducing  hqs. 
(H)— (20)  into  Eq.  (9),  whereby 

(a-  - v)t'Grr  ~ ~T  1 [‘-'re  ~ ~ ^2)] 


r 
/I 


(21) 


(a-  - ir)-^-^Gee  G^ ) - (ir  • v-)ki  (vdlj,  - 


L = 0 


and 


L = - 


4uv 


[tr  — r„  cos  6)  (Gj  — F_)  ' r r*  sm  6Gr] 


^ 2(u2  - v^) 


(22) 


|^(r  - r^  cos  6)Gj.  - sin  6(C>(,  - r.,)J 

— (Fj  F2)-^|2(u^  - v^)(r  - r^  cos  9>v„  sin  9 2uv(l  - 2r^  sin^^j 

— E,[a^f  — (r„v  sui  9 r„u  cos  9 — ur)^]f 


where 


E,  s 


^ 1 - MF 


' 1 - .Mi  sm'/3 


E, 


- sin  2^ 


(24) 


1 - Mi  sinV 

and  f and  p are  defined  in  Kqs.  (15)  and  (18),  respectively.  The  radial  and  circumrerential 
velocity  components,  v and  u become 


-I ) r/-- 


— - (rj  Fj)Ejr,,  sin  e 


1 


M 


(25i 


r)i.:,(r,.eose-r)  u,; 


(20) 


where 


'd  = y [f’"  c.-os(e 

G - ko)  - 2r  r,c  cos(q  - k-) 

I'i  C(Js(d  - G 

iO] 

(271 

I'l  - ^ [r-  sin(e 

Q - !•(;)'  2rr,.  ,sin(G  - k,)  - 

I'i  sin  (6  - G 

k,)] 

(28) 

The  boundarv  condition  of  the  vanishing  of  normal  velocity  on  the  surface  of  the  airfoils. 
Eq.  (11)  becomes 

Gp  = ^ [(Fj  Fj)  Eji',^  sin  0 - v’l]  at  r r^  and  r=] 

The  Kutta  condition.  Ecj.  (12)  (jecomes 


(2!)) 


1 


1.'. 


rjEi(r^  cos  5 — r)  F;, 


^ Ei(r„  cos  e - r) 


r 


^ ' *^'1 


(30) 


at  each  trailing  edge.  This  gives  to  two  linear  expressions  for  Fj  and  Fj  which  can  be  solvefi 
at  each  step  m the  iterative  numerical  solution. 

An  additional  Ixiundary  condition  must  be  enforced  at  the  point  representing  infinity.  The 
entire  solution  there  is  given  by  the  sum  of  the  free  stream  and  circulatory'  flows.  No  addi- 
tional velocities  must  be  generated  In’  reduced  potential  G.  So,  in  effect,  G must  be  a con- 
stant at  infinity  and,  therefore,  a constant  in  all  directions  at  r - r,,,  6 = 0.  For  conve- 
nience this  constant  is  chosen  as  zero. 

Solutions  to  the  boundary  value  problem  comprising  Fqs.  (10),  (21)  — (30),  can  be  obtamed 
by  iteration  and  the  numerical  scheme  will  Ire  rlescribed  in  the  followmg  section.  A start- 
ing point  for  the  iteration  process  is  provided  lay  Legally ’s  solution  (Ref.  13)  for  the  incom- 
pressible flow  over  two  circular  cylinders,  which  in  ternrs  of  the  reducerl  potential  G(r,  6) 
mav  be  writte.i  as 


G(r,  6>  = Real  kj 


1 r 30  Z 


i=i 


1 ^ 


> gl(t  -a) 


,-2J 


J = 1 


1 - 


r. 


(31) 


- i(ri  F^) 


- log(l  - 


) r,o  / 

(l  - r,o/'.,r'0( 

1 - r'’^) 

' *-»^G  ' 
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4.  NT'MhRICAL  l OK.Ml  LATiON 


Tile  Neumann  Ixnindaiy  value  problem  I'ornuilatetl  in  the  preeediiig  sections  is  solved 
numerically  using  the  ideas  and  techniques  developed  lor  mixed  sulisonic  and  supersonic 
riovvs  about  a single  airfoil.  The  reduced  potential  equation,  Eq.  (21),  is  solved  by  a suc- 
cessive column-relaxation  algorithm  utilir’ing  lype-dependent  dilferencing  originated  by 
Murman  and  Cole  (Ref.  2).  However,  the  sequences  ot  mappmgs  produces  a computtitional 
plane  in  which  the  coordinate  Imes  are  not  aligned  with  the  streamwise  direction.  Jameson  s 
technique  (Ref.  14)  for  developing  a cooniinate  invariant  or  “rotated”  difference  scheme 
proves  valuable  in  the  numerical  formulation  of  this  problem. 

The  solution  of  the  flow  field  is  obttiined  by  re|)lacing  the  potential  equation  with  a finite- 
difference  analog  at  each  point  in  the  grid.  The  resulting  set  ol  diflerence  equations  is 
solved  by  iteration,  a process  which  can  l^e  regarded  as  a sequence  ol  steps  in  an  artiiicial 
time  coordinate  frame.  During  each  iteration  tlie  updating  of  an  initi.-d  potential  fiehi  is  ac- 
complished simulttuieously  at  a numljer  of  |X)ints  on  the  same  column  (X  = constant)  and 
sweeping,  column  by  column,  through  the  computational  grid.  It  is  imperative  that  the  How 
field  be  swept  in  a direction  which  is  no  more  than  9tr  Irom  the  streamwise  direction  in  or- 
der to  retiiin  the  proper  domain  of  dependence  everywhere.  .A  scheme  whicli  has  pro\eil 
suitable  in  this  problem  has  been  to  divide  the  How  field  in  two  !jy  the  coordinate  ring  |)ass- 
ing  through  the  infinity  point  (r  - ly  or  V = "i.  Eacli  hall  ol  the  total  Held  can  tlien  be  con- 
sidered to  surround  a single  clement  and  is  swept  Irom  the  Icailing  edge  to  the  trailing  edge 
of  the  airfoil  along  one  surface,  and  tlien  along  the  other  as  illustrated  in  Tigs.  2 and  3. 
column  of  points  extending  from  the  surface  to  the  middle  coordiiuite  ring  is  thus  updatial  at 
each  pass.  .At  the  end  of  each  sweep  the  Kutta  eomlition  at  each  airfoil  trailing  edge  is  en- 
forced  by  solving  for  the  circuhition  values.  E,  and  Fi.  which  salisly  tlie  diffi'rence  equa- 
tions corresponding  to  Eq.  (3(f). 

Hefore  Eq.  (2D  can  be  used  for  numerical  simulation  it  shouh!  oe  re-wi  itlen  in  terms 
of  the  working  coordinates  X and  A’.  .Since  6 0(X)  and  r r(A')  flcfined  in  lAis.  t(3)  — (S). 

(•  - \-'c 

(■9  - (-X  - ‘'X 


where  i and  j are  the  indices  in  the  X and  Y riirections,  respectively,  and  the  superscript 
denotes  “old”  values.  A similar  expression  can  he  wi'iltcn  lor  G^,.  As  a result,  velocities 
are  frozen  at  values  of  the  [)otential  based  on  the  previous  sweep.  Second-order  derivative 
terms  are  written  takinj^  into  account  their  domain  of  dependence.  .\t  sulrsonic  points  cen- 
tral differences  are  appropriate.  The  latest  values  of  the  potential  are  used  with  the  cur- 
rent column  (j)  beiiifi  relaxed.  Thus,  if  one  denotes  the  difference  between  the  new  :ind  old 


1 


r 


values  by  6Gj,  then  the  secoad-order  accurate  dUTerenc:e  formulas  are  written  as 


AY 


‘■XX  “ 

AX2 

■ ‘^1-1, j 

-(0?, 

JO  - 

^^1,  JO 

- 26Gj^j  <^^’j,j-j) 

( k 


“ 4AXAY 
where  oi  is  the  relaxation  factor. 


I ^'1-1,  j+i 


^"'1-1,1-!) 


At  supersonic  points  the  proper  domain  of  dependence  extends  upstream.  Following 
Jameson,  Ref.  14,  the  second  derivative  terms  are  broken  up  into  their  streamwise  and 
normal  components  as 

^^2  2 ,^2 

^SS  ■ ~7l  ^^2^’XY  ■ ^3*^'YY 


and 


where 


*-■83  “ i^l^’XX  * i^2^'XY  ' I^3*-'YY 


G 


'NN  ~ ^'lt'.xX 


B, 


CRixY  ^'^‘-'YY 


irfX'- 


B)  = 2uvf 


X'Y' 


(.>3) 


IG  = vfV'- 


t'l  V f 


Co  — 3uv 


X'' 

1' 

r 


G^  = irfV'' 

The  derivatives  mal<ing  up  the  normal  eomi'onent.  Gxn-  approximated  by  the  c-entral 
difference  formulas  given  above  with  the  rela_\ation  factor  taken  as  otu'.  l ])wind  differencing 
is  used  in  both  the  X and  Y directions  to  evaluate  tiu-  sc'veral  parts  of  Gs-s,  the  streamw  ise 
comptinent.  Thus,  for  examide,  if  u • n and  v i)  the  contribulitms  will  !)(■ 
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^•xx 


AX 


G.?,,_  > 26Gj) 


1!) 


26C5j  - 26Gj.j) 
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YY 


AY 


,2 


(G°  , - 2G°.,., 


no 
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" 2AXAY 


(G°i, , - Gi-i, , 


26Gj  - 6Gj-i) 


Similar  expressions  can  be  written  for  cases  where  u and  v take  on  different  signs.  It  should 
be  noted  that  Jameson’s  rules  for  balancing  “old”  and  “new”  values  to  provide  the  right  arti- 
ficial viscosity  have  been  followed. 


In  order  to  satisfy  the  boiuidary  condition  at  each  of  the  airfoil  surfaces  a dummy  row  of 
points  is  created  beyond  each  surface.  The  requiremeiit  that  there  be  no  flow  through  the 
surface  at  each  of  the  surface  grid  points  (j  = jj)  can  then  be  met  by  requiring  that  the 
change  in  potential  at  the  dummy  point  equal  the  change  at  the  first  point  off  the  surface, 
that  is,  6Gj  = hG,  .j.  By  this  device  4>r  at  each  surface  will  remain  zero.  At  the  mid-ring 
boundary,  which  has  been  set  up  to  divide  the  various  sweep  regions,  the  values  of  the  po- 
tential from  the  previous  pass  are  kept;  essentially  this  is  a Dirichlet-type  boundarj’  condi- 
tion. 


By  applying  these  finite-difference  relations  simultaneously  at  all  grid  points  on  a 
column  extending  from  either  airfoil  surface  to  the  mid-ring,  a system  of  linear  algebraic 
equations  for  the  corrections,  6Gj  is  obtained.  This  system  can  be  represented  by 

[Pil][6Gj]=  [Kj 

It  can  be  easily  checked  that  the  matrix  [Pjj]  is  tridiagonal  and,  therefore,  can  be  in- 
verted easily.  Sequentially,  column  by  column,  all  |X)ints  in  the  fields  are  updated,  except 
for  the  point  representing  infinity  which  is  kept  at  its  initial  value. 

At  the  end  of  each  sweep  new  values  of  the  circulation  can  be  computed  by  applying  Eq. 
(30)  at  each  trailing  edge  and  solving  the  two  resulting  equations  fur  Tj  and  To.  But  this  ex- 
pression, which  can  be  rewritten  concisely  as  <l>9  = U,  was  found  to  be  the  source  of  errors 
in  some  cases.  Because  of  the  complicated  geometries  being  dealt  with  in  this  problem, 
large  circumferential  gradients  of  the  potential  could  appear  at  one  of  the  trailing  edges. 

For  example,  there  could  be  a shock  or  an  expansion  fan,  or  both,  veiT  near  the  trailing 
edge  of  a slat.  For  these  cases  we  find  the  stability  to  lie  enhanced  by  under  relaxation  of  the 
evaluation  of  the  circulation  constants. 
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The  differencing  pi’ocedure  described  above  is  of  a noni'onservative  form.  Although  this 
method  produces  inaccurate  shock  “jump”  conditions  (Reis.  15,  16i  generally  Ijctter  agree- 
ment with  experimental  data  is  obtained.  This  is  likely  lortuitous.  A completely  satislac- 
torj-  solution  would  be  obtained  only  with  a fully  conservative  treatment  of  the  inviscid  equa- 
tions together  with  a proper  modeling  of  the  viscous  effects  and  of  tlie  shock  wave-boundaty 
layer  interaction  process  in  particular. 

To  increase  diagonal  dominance  in  the  tridiagonal  matrix,  additional  artificial  viscosity 
can  be  added  in  the  form  of  a term  tG„t>  t is  the  lime-like  coordinate  (Kef.  14).  With 

the  basic  tinie-like  step  being  the  iteration  cycle  this  additional  term  can  be  represcnteti  by 

A/  / 1 II  , V \ 


il  X.,-  1 


For  V ■ 0 this  can  be  put  in  the  form 

\ r AX  \ * AX  \ ' ,,  r//'  i 'O  \i^ 

^ AY  X' ~ L '' AY  ~ ‘^'1-1.  j)j^' 

This  term  is  added  to  the  finite-difference  analog  of  Hq.  (lil).  As  discussed  in  Ucl.  14,  e 
should  be  only  large  enough  to  ensure  convergence.  In  practice,  e has  Ijccn  taken  to  lie  zero 
in  most  cases. 

Finally,  in  order  to  improve  the  computational  efficiency,  the  calculations  are  per- 
formed on  a sequence  of  two  grids.  The  initial  (coarse)  grid  is  computcil  starting  with  an 
initial  representation  of  the  reduced  potential  G(X,  Y)  given  Ijy  Fq.  (31).  The  starting  con- 
ditions for  the  second  (fine)  grid  is  obtained  by  a linear  interpolation  of  the  converged  coarse 
grid  solution. 
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5.  UESll.TS 

The  method  has  been  applied  to  a variety  of  two-element  airfoil  configurations.  Figures 
(>,  7,  8 show  the  c(miputed  flow  field  about  a NACA  23012  airfoil  with  a 2-11  flap  at  M =0.5 
and  an  angle  of  attack  of  5°.  The  calculation  was  performed  using  a series  of  two  grids  with 
the  final  one  being  00  (in  the  circumferential  direction)  by  30  (in  the  radial  direction).  The 
computed  pressure  distriljutions  on  the  main  element  and  flap  are  shown  in  Fig.  6,  and  the 
computed  streamline  patterns  in  the  physical  and  annular  domain  arc  shown  in  Figs.  7 and  8. 
Imbedded  regions  of  supersonic  flow  can  be  seen  on  the  latter  graphs.  In  Fig.  0,  the  large, 
“tic"  mark  on  the  Cp  axis  denotes  sonic  pressure  and  the  sharp  nose  sliock  computed  is  evi- 
dent. 

In  this  calculation  the  changes  in  the  reduced  potential  were  converged  to  3 10"’’  on  the 

final  grid  in  approximately  five  minutes  on  an  IBM  370/168  computer.  Some  preliminary 
studies  have  indicated  that  the  application  of  the  eigenvalue  extrapolation  method,  as  dis- 
cussed in  Ref.  8,  could  substantially  improve  convergence. 

The  solution  for  the  flow  over  a Clark  Y airfoil  with  a 30  ( Miixwcll  slat  with  a gap  height 
of  10  0 at  M*  = 0.6  is  shown  in  Fig.  9.  The  supersonic  regions  in  this  case  arc  considerably 
larger.  On  the  slat  the  entire  upper  surface  flow  is  supersonic.  This  results  in  a large  lift 
being  carried  by  the  slat  with  a significant  amount  of  "rear"  loading.  The  streamline  and 
sonic  line  patterns  in  the  physical  and  annular  domains  are  sliown  in  Figs.  10  and  11. 

A still  larger  supersonic  region  occurs  in  the  next  example  which  shows  a slatted  con- 
figuration that  is  olitained  from  a modified  N.\CA  (i-lAUdS'  section  profile  at  M,  0.7  and 
a = 6°.  llie  computed  jjressure  distribution  and  streamline  contours  aia  sliown  in  Figs. 

12  and  13,  Even  though  the  present  grid  is  not  as  densi  as  thosi  used  in  Mu  single-element 
computations,  the  shock  is  clearly  defined.  It  should  Ik'  kept  in  mind  that  the  fkiw  on  the 
lower  surface  of  the  slat  could  separatt  so  that  sonu  of  tlu  lift  shown  Ik  ing  carried  by  the 
slat  could  be  lost.  A calculation  for  the  M„  o.",  J flow  ovi  r a slatUd  configuration 
based  on  a modified  NACA  t>4A009  airfoil  is  piusi  iUetl  in  Fig.  11.  A plot  of  the  Mach  num- 
ber contours  illustrating  the  large  regions  of  supersonic  flow  is  given  in  Fig.  15. 

•)0 
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Fig.  6 Computed  Surface  Pressure  Distribution  NACA  230*12  Airfoil/2H  Flap,  M 


MAIN  AIRFOl  L 


Fig.  7 Computed  Streamlines  — Annular  Domain  NACA  23012  Airfoil/2H  Flap,  Moc  = 0.5,  0 = 4° 

The  choice  of  sweep  directions  causes  no  difficulties  for  cases  where  a supersonic  re- 
gion should  span  the  entire  gap  between  the  two  elements.  This  situation  occurs  for  the 
slatted  NACA  0012  airfoil  at  =0.6  and  a = 4°  as  seen  in  Fig.  16.  Here  Mach  number 
contours  are  plotted  with  the  sonic  line  being  marked  more  heavily.  The  corresponding 
pressures  are  shown  in  Fig.  17.  The  supersonic  region  spanning  the  gap  is  swept  in  the 
streamline  direction  along  both  surfaces  so  that  the  proper  domain  of  dependence  is  main- 
tained. 

A final  example  provides  some  measure  of  the  performance  of  the  method.  Figure  18 
depicts  the  computed  pressures  and  experimental  data  for  the  geometry  shown.  The  geome- 


PHYSICAL  PLANE 


Fig.  8 Computed  Streamlines  — Physical  Plane  NACA  23012  Airfoil/2H  Flap,  ~ 0.5,  o - 4 


try  used  in  tlic  computation  is  the  (i4A40()  airfoil  witli  7.  8 A slat  configuration  testtd.  cor- 
rected through  the  addition  of  a nominal  displacement  thickness.  The  overall  agreement  is 
not  very  good  and  this  implies  that  the  Ijoundarj’  layer  correction  applied  to  the  ordinates  was 
insufficient.  The  excellent  agreement  and  the  lower  surface  suggests  that  discrepancies  on 
the  upper  surface  of  the  main  airfoil  and  on  the  slat  are  due  to  the  trailing-edge  separation 
on  the  main  foil  and  separation  on  the  lower  surface  of  the  slat.  The  disagreement  between 
the  numerical  solution  and  the  data  on  the  upper  surface  of  the  main  element  consists  mainly 
of  an  upward  shift.  Despite'  the  shift  all  thefe'atures  of  tlie  flow  are  predicted  by  the  computa- 
tion including  the  double  peak  occurring  near  the  nose  of  the  main  foil.  The  agreement  would 
undoubtedly  be  improved  by  profierly  accounting  for  viscous  effects.  .-Xs  shown  in  Fig.  19, 
the  sonic  region  not  only  spans  the  entire  gap  but  also  touches  the  trailing  edge  of  the  slat. 
With  the  particular  differencing  scheme  used  in  the  application  of  the  Kutta  conditions  the 
method  had  no  difficulty  computing  the  rapid  e'X])ansion  around  the  trailing  edge  of  the  slat. 


Fig.  9 Computed  Surface  Pressure  Distribution  Clark  Y Airfoil/30%  Maxwell  Slat,  M, 


MAIN  AIRFOIL  X 
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Fig.  10  Computed  Streamlines  — Annular  Domain  Clark  Y Airfoil/30'/o  Maxwell  Slat,  M 


= 0.6,  a = 6' 


Fig.  11  Computed  Streamlines  — Physical  Plane  Clark  Y Airfoil/30%  Maxwell  Slat,  M . = 0.6,  Q = 6 


4 + + 4 + 


MAIN  AIRFOIL 
SLAT 


Fig.  14  Computed  Surface  Pressure  Distribution  Modified  NAC.  A009  Airfoil/Slat.  M^  = 0.8,  q = 4 
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Fig.  15  Mach  Number  Contours  Modified  NACA  64A009  Airfoil/Slat. 


Fig  16  Mach  Number  Contours  NACA  0012  Airfoil  and  Slat, 


Fig.  18  Computed  Surface  Pressure  Distribution  NACA  64A406/7.8A  Slat,  M - 0.649, 


G.  CONCLUSIONS 
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A numerical  method  has  been  develoj^ed  to  calculate  the  inviscid  transonic  flow  over  ar- 
bitrary two-clement  airfoil  systems.  Conformal  mapping'  procedures  were  utilized  to  de- 
velop a convenient  computational  domain  where  boundary  conditions  on  the  airfoil  surfaces 
and  at  infinity  could  be  accurately  prescribed  and  witli  appropriate  spatial  grid  definition  to 
account  for  steep  flow  gradients,  'llie  exact  irrotational  flow  equations  were  derived  in  terms 
of  a smoothl}'  varjing  single-valued  reduced  potential  in  the  computational  domain.  A stable, 
accurate  and  efficient  finite-difference  procedure  was  developed  using  mixed-flow  relaxation 
techniques.  The  numerical  method  was  applied  to  a variety  of  airfoils  with  leading-edge 
slats  and  trailing -edge  flaps  over  a range  of  free-stream  lUacli  numbers  and  angles  of  at- 
tack. 

Our  computed  results  indicate'  the  ability  of  the  method  to  compute  inviscid  transonic 
flows  over  two-clement  airfoil  configurations.  One  aspect  of  this  effort  which  requires  fur- 
ther study  is  the  improvement  of  the  computational  efficiency.  We  have  preliminary  com- 
putations with  tlie  eigenvalue  extrapolation  procedure  of  Ref.  8,  which  indicate-  a possible 
40  0 reduction  of  iteration  cycles.  'Tliis  technique  will  be  implemented  in  the  near  future. 

Another  area  of  further  study  is  the  verification  of  the  accuracy  of  our  method  by  com- 
parison with  experimental  data.  However,  it  is  evident  tiiat  transonic  two-element  airfoil 
flows  arc  strongly  influenced  by  viscous  effects.  Thus,  we  are  in  the  process  of  including 
boundary  layer  effects  in  our  computational  techniciue.  Special  emphasis  will  Ije  given  to  in- 
cluding strong  viscous  interaction  effects  at  trailing  edges,  shock  waves,  and  in  regions  of 
separated  flow. 

We  also  intend  to  reformulate  tlic  proljlem  in  the  design  (inverse)  mode,  wliereliy  the 
pressure  distribution  will  be  prescribed  over  all  or  parts  of  tlie  airfoil  surfaces  and  the  or- 
dinates will  evolve  as  part  of  tlie  solution  procedure.  Technicjues  devclopt'd  for  single-ele- 
ment design,  sucli  as  Ref.  17,  will  be  utilized. 
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APPENDIX 

Evaluation  of  the  Z^— Z4  Mapping  Coefficients 


VVe  discuss  in  detail  the  mapping  from  the  Z3  to  tlie  Z4  phme  as  given  in  Eq.  (3).  In  or- 
der to  determine  tlie  coefficients  of  this  mapping  we  consider  the  following  preliminary  map- 
ping 

g =f'/^ 


(A.  la) 


g = 


f = 


(A.  lb) 


(z  - b)(z  - i) 
where  c and  d are  the  unknown  coefficients. 

On  the  unit  circle  zz  = 1,  f may  be  written  as 

f 


(A.  Ic) 


z - a 


I2 

a 


z - b 

so  that  the  unit  circle  in  the  z plane  maps  to  a slit  in  tlte  f plane  at  an  angle  arg  (b/a)  with 

z - a 


end  points  at  fj  and  lo  where 


z — b 


reaches  its  extreme  values.  Then  through  the  map- 
ping (A,  la)  we  get  another  slit  stretching  from  (f))’  '*  to  (f,)'  Now  the  mapping  (A.  lb) 
maps  a unit  circle  in  the  < plane  (i'^  = l)  to  a slit  Ijetwcen  gq  and  gj.  Thus,  if  we  choose  c 
and  d so  that 


<r  f ^ ^ 


1 


(A.  2) 


g2  =12'" 

then  the  unit  circle  in  the  z plane  will  be  mapped  to  a unit  circU'  in  the  c plane.  To  make 
the  slits  collincar  requii’cs 

a.-g(2)  f (A.3) 

Then  c and  d must  be  chosen  so  that  g,  and  g.,  have  tlie  magnitude  appropriate  to  E(|.  (A.  2). 

To  determine  g,  and  go  take  c 1/c,  d = l/d  and  evaluate 
]_  dg  1 ^ 1 1 1 c-d  c-ci 


g (If 


f - c 


- c 


r - d 


- d 


(f-e)  ('  -d)  (i-c)  (f-d) 


V = 0 


(A. 4) 


,'i.S 


J 


The  ends  of  the  slitiu'ctlie  brnnch  points  gj  nnd  g2  where  vanishes.  Thus,  at  a branch 
point  £ satisfies 


0 = (e  - cKc  - d)(c  - d)  - i?  - c)(f  - d)(c  - d) 
= - Bt  + C 


( A . 5 a) 


where 


A - (c  - c)  - (d  -■  d)  (A.  5b) 

1)  2(cc  - del)  (A.  5c) 

C = cc{d  4 cl)  — del  (c  ‘ c)  (A.5d) 

Then  if  f,  and  t,.,  are  the  roots  of  the  branch  points,  g,  and  go  arc  given  iiy 

ft  ■ 

The  cxiJrt'ssions  for  tiec  Ijranclt  points  fj,  fo  may  be  ewaluatc'd  in  an  anaiagous  manner.  Now 


from  condition  (A.  la) 


= (I'lf;;)'' 


Ttilizing  Eets.  (A. 5)  we  obtain  ' jeo  - ^ mid  £i  + tA  " \ with  Kei.  (A.  (>)  we  can  develop 

A 

(after  some  tedious  algelira)  tliat 

/c  - cV  RlcI-  - l)  d1'  .. 

Vd  - el  ^ [(I  dl  " - l)  cj 

and  an  anaiagous  expression  few  f, f . Thus.  Ih).  (A. 7)  rceiuircs  that 

- ~ l'  (A.  9) 

d - d VI)  - b I 

Now  till’  lirancli  point  condition  l.ii.  (A. 4)  ma\'  be  rewritten  as 

t2?  - c - e')  (d<“  - d - d)  , . 

o — ; , ;-7 ^ (A.  10) 

(t  - c)(r  - c)  (?  - eO(f  - d) 

and  the  branch  poi.it  gi  given  in  Ki|.  (A.ii)  lucomes 

ft 

However  £j  lies  on  tlu'  unit  circle,  and  w i are  tree  to  rot.ite  the  ; iilane  so  that  the  ai'g  (£j) 
talces  any  desired  value.  Thus,  if  we  assign  .i  value  for  < ] . tlu'  condition 

g,  (f,)‘  ' (A.  12) 

givt'S  one  relationshi()  foi'  c,  d.  c,  d .md  l.i|.  (.\.  9)  give's  a se'cond;  :uid  since*  c and  d are  de*' 
termined  from  c and  d we  ha\e'  two  e'eiuations  for  c and  d. 
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These  equations  may  be  developed  in  a convenient  form  by  noting  that 


f, 


= (zi  - a)(zi  - a)  ^ 


‘ (zi  - b)(zi  - b) 
while  according  to  (A.  9)  we  define  a quantity  E as 


z,  — a 

Zi  - b 


2 b 


^ c - c _ /a  - aV''*'  _ a 1 2 - i \ b 
^ " d-d  "M)  - b/  ~ /\|b|2  - 1 / a( 


(A. 13) 


{A.  14) 


Thus 


Z)  — a j 2 ' i !j  1 2 — 1 
Zi  - b V|  al  2 - l/(  ^ 


= PE 


(A.  15) 


The  two  equations  for  c and  d are  obtained  from  Eqs.  (A.  9)  and  (A.  10)  utilizing  Etjs.  (A.  14) 
and  (A.  15)  as 


-^44) 


- j) 


(A.  Id) 
(A.  17) 


where  D 2fj(l  — pE)  and  £,  has  an  assigned  value  on  the  unit  circle.  Combining  the  above 
two  equations  we  obtain 


2c  Dt^E(,l+i)*E(d-i) 
| = D*dE(d*i)-.:(cl-i) 


(A.  IM 
(A.  19) 


We  can  also  note  that  i)y  definition 


/ 

1 \ 

arg(c) 

argl^c 

;+rj  = arg(c 

' “c) 

arg  (c 

l+^)  = arg(c 

'4) 

1 = arg(d) 

wherel)y  E(|S.  (A.  10)  and  (A.  17)  yield 

arg(c)  ~ arg(Ed)  arg(Dl  (A.  20) 

Taking  the  complex  conjugate  of  E(i.  (A.  18)  and  multiplying  by  Fai.  (A.  19)  and  simplifying 
the  exjn'ession  using  Ecj.  (A.  20)  gives  the  folloyving 


Letting  r 


I I)  I I E I 


d + yve  obtain 
d 


I E I 


^'-d 


4 


(11)1  . r)2  - r-  ‘ 4 IE  I ■ 1 


(A.  21) 


(A.  22) 


10 


which  may  easily  be  solvecl  since  1)  and  K are  known.  Ihen 

1 r 


Idl  + 


d 


1 E 


wheretjy 


Idl  = 


r" 


TTI  - 1 


(A.  23) 


2 1 E I V 4 1 E 

where  the  root  larger  than  1 has  been  cliosen.  Also  from  Eq,  (A. 20) 

ai'g(d)  - arg(D/E)  (A. 24) 

and  thus  Eqs.  (A.  23)  and  (A.  24)  determine  d and  c may  then  1)c  found  from  Eq.  (A.  10). 

The  mapping  coefficients  found  completely  determine  the  mapping  given  by  Eqs.  (A.  1). 
This  mapping  differs  from  the  one  developed  l)y  Ives  (Kef.  11)  and  given  in  Eq.  (3)  by  a scale 
transformation.  If  we  define 


S K 


■ 1 I ^ 

c + - - d - 
c ^ 

1 , i 

a + - - l3  - I 
a b -1 


(A. 251 


Then  Eqs.  (A.  1)  arc  identical  to  Eq.  (3)  with 

E.i  = £/S 
^■1T  ~ c/S 

Zjj;  = d/s 


and 


1-  = SS 


Z3  - /. 


/,  ,1  Cl 

K o ~ K 


Thus  we  have  reduced  tlie  determination  of  the  mapping  coefficients  of  Ec(.  (3)  to  the  solution 
of  two  real  (luadrutic  eciuations.  In  Kef.  11.  these  mapping  coefficients  are  found  using  a 
complicated  three-dimensional  Ncwton-Haiiheson  iteration  lu'oce.ss. 
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•Vir  X'ehicle  Technology’  Dept. 

ATTN:  Code  301  1 


Director,  Ames  Directorate 

F.  S.  .Vrmy  .\ir  Mobility  RitD  Lab. 

■Vines  Research  Center 

Moffett  Field,  CA  94035  1 
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Director,  Langley  Directorate 

Lo  So  Army  Air  Mobility'  R&D  Lab. 

Langley  Research  Center 

Hampton,  VA  23665  1 

Director,  Eustis  Directorate 

U.  So  Army  Air  Mobility,’  R&D  Lab. 

Fort  Eustis,  VA  23604  1 

V.  So  Air  Force  Flight  Dynamics 
Laboratory 

Wright- Patters  on  AFB,  OH  45433 
ATTN:  PR,  Prototype  Division  1 

FXM,  Aeromechanics  Branch  1 

Air  Force  Office  of  Scientific  Research 
Directorate  of  Aerospace  Sciences 
Bolling  AFB,  DC  20332  1 

National  Aeronautics  and  Space 
Administration 
600  Independence  Avenue,  SW 
Washington,  DC  20546 
ATTN:  Code  R^VA  1 

Code  RAV  1 


Lockheed  Missiles  & Space  Co. , Inc. 
Huntsville  Research  & Engineering 
Center 

P.  O,  Box  1103 
Huntsville,  AL  35807 

ATTN:  Mr.  A.  Zalay  1 

Rockwell  InterTiational 

Science  Center 

Thousaird  Oaks,  CA  91360 

ATTN:  Dr.  N.  Malmuth  1 

Nielsen  Errgiireering  & Research,  Inc.  1 

510  Clyde  Avenue 

Mountain  View,  CA  94043  1 

University  of  Cincinnati 
Department  of  Aerospace  Engineeriirg 
and  Applied  Mechanics 
Cincinnati,  OH  45221 
ATTN:  Dr.  R.  T.  Davis  1 


National  Aeronautics  and  Space 
Administration 
Ames  Research  C''nter 
Moffett  Field,  CA  94035 
ATTN:  Dr.  T.  Gregory  1 

Dr.  G.  Chapman  1 

National  Aerodynamics  and  Space 
Administration 
Langley  Research  Center 
Hampton,  VA  23665 

Subsonic,  Transonic  Aerodynamics  Div. 
ATTN:  Dr.  lames  F.  Campbell  1 

Dr.  Stephen  Womom  1 


Flow  Research,  Inc. 

1819  South  Central  Avenue 
Kent,  WA  98031 

ATTN;  Dr.  F.  Murman  1 
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